Flow diagram of the metal-insulator transition in two dimensions 
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The discovery of the metal-insulator transi- 
tion (MIT) in two-dimensional (2D) electron sys- 
tems [1] challenged the veracity of one of the most 
influential conjectures [2] in the physics of disor- 
dered electrons, which states that "in two dimen- 
sions, there is no true metallic behaviour"; no 
matter how weak the disorder, electrons would be 
trapped and unable to conduct a current. How- 
ever, that theory did not account for interactions 
between the electrons. Here we investigate the in- 
terplay between the electron-electron interactions 
and disorder near the MIT using simultaneous 
measurements of electrical resistivity and mag- 
netoconductance. We show that both the resis- 
tance and interaction amplitude exhibit a fan-like 
spread as the MIT is crossed. From these data we 
construct a resistance-interaction flow diagram of 
the MIT that clearly reveals a quantum critical 
point, as predicted by the two-parameter scaling 
theory [3]. The metallic side of this diagram is ac- 
curately described by the renormalization group 
theory [4] without any fitting parameters. In par- 
ticular, the metallic temperature dependence of 
the resistance sets in when the interaction ampli- 
tude reaches 72 w 0.45 — a value in remarkable 
agreement with the one predicted by theory Q]. 

The low amount of disorder in high mobility sili- 
con metal-oxide-semiconductor field-effect transistors (Si 
MOSFETs) allows measurements to be made in the 
regime of very low electron densities where correlation 
effects due to electron-electron interactions become es- 
pecially important. (Ratios rg = Ec/Ep > 10 be- 
tween Coulomb and Fermi energies are easily reached 
with Fermi energies of the order 0.7 meV.) This mate- 
rial system has the additional advantage that its elec- 
tron spectrum has two almost degenerate valleys, which 
further enhances the correlation effects. Indeed, the low- 
temperature drop of the resistance on the metallic side 
of the transition (FigHlJa)) in Si MOSFETs is the most 
pronounced among all 2D electron systems [1, 0] . 

At low temperatures fc^T < h/r < Ep, electrons prop- 
agate diffusively. Here, T is temperature, r is the electron 
mean free time, and Ep is the Fermi energy. This region, 
referred to as the diffusive regime, extends up to a few 
kelvin near the MIT. As a result, the low temperature 
physics of the disordered electron liquid, combined with 
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FIG. 1: Temperature dependences of the resistivity and of 
the strength of spin-related interactions for electron densities 
across the MIT. (a): p{T) traces at B = 0. (b): Cee{T), ex- 
tracted from the slope of the magnetoconductance (see text). 
Both p{T) and Cee{T) display a fan-like structure across the 
MIT. The densities (in units of 10^° cm"^) are: 7.24, 7.53, 
7.83, 8.26, 8.70, 9.14, 9.58, and 9.87. In panel (a), the den- 
sities increase from top to bottom, while in panel (b), they 
increase from bottom to top, i.e., the interaction increases 
in the metallic phase and decreases in the insulating phase as 
the temperature is lowered. M and I regions indicate metallic 
and insulating phases, respectively. 



the large values near the MIT, is determined by the 
properties of interacting, diffusing electrons 0, B Q- 

With this in mind, and with the view of studying 
the role of disorder and effective strengths of the elec- 
tron interactions in the vicinity of the MIT, we have 
performed simultaneous measurements of the resistiv- 
ity, p, and in-plane magnetoconductance, cr{B). Since 
the origin of the magnetoconductance lies in the inter- 
action corrections to the conductance involving differ- 
ent spin projections (Hartree-like) , the corresponding in- 
teraction strength can be determined from the slope of 
(j{B). When a magnetic field B is_ applied parallel to the 
2D plane in the diffusive regime [loillli, the magnetocon- 
ductance Aa{B,T) = a{B,T) - a{0,T) is proportional 
to in the weak field limit b = g^sB /ksT <C 1. On the 
other hand, at higher temperatures K/t < ksT < Ep, 
referred to as the ballistic regime, Aa{B,T) is expected 
to be proportional to Tb^ ^] . Because the magnetocon- 
ductivities in these two regions are so different, measuring 
Aa{B, T) provides a reliable way to identify the diffusive 
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region. 

The samples we used had peak electron mobilities 
about 3 X 10* cm^/Vs at 0.2 K. The resistivity p was 
measured by a standard low-frequency lock-in technique. 
The magnetoresistance traces, p{B)^ for a representative 
density in the diffusive region rig = 9.14 x 10^° cm~^, are 
plotted in Fig. [H^a) for different temperatures; p{T) at 
B = for this density is the third curve from the bottom 
in Fig. [TJa). The corresponding magnetoconductivities, 
Act — 1/ p{B) — l/p(0), are plotted as function of &^ in 
Fig. IHb). When plotted this way, one can see that the 
(t(6^) curves are linear. It can also be seen that the slopes 
decrease slowly with temperature. This is important, be- 
cause had the electrons been in the ballistic regime, the 
slopes would have instead increased with temperature by 
an order of magnitude in the temperature interval used 
due to the expected dependence. By applying the 
same procedure to different electron densities, we find 
that the diffusive region extends to approximately 25% 
in electron densities above the critical density of the MIT. 
(The critical density rtc ~ 8 x 10^" cm^'^ in this sample.) 

The explicit form of the slope of A(t(6^) for the single 
valley case was derived in Refs. 111] . It is straight- 
forward to accommodate the valley degrees of freedom. 
In the case when there are degenerate valleys (n„ = 2 
for silicon), we obtain for Act (6^), in the limit 6^1, the 
expression 
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Act(62) = -0.091— nf, 72 (72 -f 1) 6^ (1) 
■nh 

(We note that the constant 0.091 in Eq. ([T]) provides a 
more accurate coefficient compared to the value 0.084 
given in Ref. 10].) Here, 72 is the effective electron- 
electron interaction amplitude in the spin-density chan- 
nel. (In standard Fermi-liquid notation, 72 is related to 
the parameter Fg as 72 = — Fg /(I -t- Fg ).) As was ob- 
served earlier in Ref. [J], the large factor due to the 
valleys enhances the effect of the e-e interactions. (In 
this paper we have assumed that the valleys are fully de- 
generate at all temperatures by restricting ourselves to 
T > T„, where T„ is the temperature scale associated 
with the valley splitting. This temperature is not pre- 
cisely known, but our analysis of the low temperature 
data suggests T„ w 0.5 K. Further details of this region 
will be presented elsewhere.) 

Eq. ([1]), as was shown in Ref. [Tlj, incorporates renor- 
malization group (RG) corrections to first order in p 
and is therefore strictly valid only deep in the metal- 
lic region. For general p, the magnetoconductance in 
the diffusive regime will still retain the 5^-form, Act — 
— {e^/irh) Cee{l2, p) n1 b^. The coefficient Cee reflects the 
strength of spin-related interactions at any value of the 
resistance (as long as gp-sB < ksT < h/r). This is be- 
cause the in-plane magnetoconductance is a consequence 
of the splitting of the spin subbands. The fluctuations in 
spin-density lead to finite temperature corrections to the 



resistivity via the electron-electron interaction amplitude 
in the spin-density channel, 72. Thus, the spin-splitting, 
by reducing these fluctuations, leads to a finite magneto- 
conductance through 72. 

The temperature dependences of the parameter Cee, 
extracted by fitting the ct(6^) traces in Fig. [U^b), are 
shown in Fig. \V[h) for various densities across the MIT. 
To the best of our knowledge, this is the first observa- 
tion of the temperature dependence of the strength of 
the electron-electron interactions. In Fig. [T{a), we plot 
p{T) at zero magnetic field for the same densities. (In 
the following, p is always expressed in units of nh/e'^ 
0.) Fig. [U reveals that not only p but also the inter- 
action strength exhibits a fan-like spread as the MIT is 
crossed. We see that, while the interaction grows in the 
metallic phase as the temperature is reduced, it is sup- 
pressed in the insulating phase. The magnitudes of these 
changes depend on how far the system is from the MIT, 
with both p and 72 becoming practically temperature in- 
dependent for densities close to the MIT. This behaviour 
is indicative of the flow around a QCP, which is fully 
consistent with the theoretical prediction that the evolu- 
tion of p and the interaction amplitude 72 with tempera- 
ture are described by a two-parameter system of coupled 
RG equations manifesting an unstable fixed point cor- 
responding to the QCP [3j]. Note that a phase diagram 
can be presented in different coordinates obtained as a 
result of various (rather broad) transformations, preserv- 
ing the topology of the phase diagram. In our case, the 
transformation involves using the slope of the magneto- 
conductance. 

To see this in a different way, we combine Figs. [UJa) 
and (b) and construct a two-parameter "fiow" diagram 
in the disorder-interaction (p - Cee) space. The flow 
diagram is presented in Fig. [31 One can see that such a 
flow diagram constructed out of p and the parameter Cee 
(which is sensitive to the amplitude 72) clearly indicates 
the existence of a QCP from which three separatrices, 
one attractive and two repulsive, can be drawn. 

On the metallic side of the flow diagram, a separatrix 
which forms the envelope of a family of non-monotonic 
curves can be clearly seen in Fig. [3] flowing toward the re- 
gion with p ^ 1. The theoretical details of the evolution 
of p and 72 in this region with p ^ 1 were discussed in 
detail in Ref. 0] in terms of an RG theory. We recall the 
salient features of the theory here. The theory predicts 
that: 

(i) The amplitude 72 increases monotonically as the tem- 
perature is reduced. 

(ii) The resistance, on the other hand, has a character- 
istic non-monotonic form changing from insulating be- 
haviour {dp/dT < 0) at high temperatures to metallic 
behaviour [dp/dT > 0) at low temperatures, with the 
maximum value pmax occurring at a crossover tempera- 
ture T = Tmax- 

(iii) The value of the amplitude 72 at T = Tmax is uni- 
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FIG. 2: Magnetoresistivity and magnetoconductivity for electron density Us — 9.14 x 10^" cm^^ measured at different 
temperatures, (a) Resistivity plotted against magnetic field applied parallel to the 2D plane, (b) Magnetoconductivity 
Aa{B,T) = a{B,T) — a{0,T) (in units of the quantum conductance) plotted as a function of (gfiBB/kBT)^ ■ The curves in 
(b) are vertically shifted for clarity. Note that the curves for different temperatures are linear when plotted vs. (gfiBB/kBT)^ 
and their slopes slowly decrease with temperature. 
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FIG. 3: The disorder-interaction (p - Cee) flow diagram of the 
2D electron liquid in Si-MOSFETs. The red dot at p ^ 0.8 
and Cee ~ 0.035 indicates the location of the QCP from which 
the three separatrices (black lines) emanate. Arrows indicate 
the direction of the flow as the temperature is lowered. Elec- 
tron densities are indicated in units of 10^'^ cm~^. 

versal in the limit pmax ^ 1, depending only on n„; for 
iiy = 1, it is 2.08, while for = 2, it has the much lower 
value 0.45. 

(iv) Although the values of pmax and Tmax are not uni- 
versal and depend on the system, the behaviours of 
p{T) j pmax and 72 (T) are universal when plotted as func- 
tions of /9maxln(T/Tmax)- 

The results of comparison between experiment and 
theory are presented in Figs. UJa) and (b), for the 



three largest densities in the diffusive regime with 
Pmax ^ 0.4 [-Kh/e^). At these densities, the resis- 
tance is small enough, and one can safely extract the 
interaction amplitudes 72 from Aa{b'^) using Eq. ([T]), 
i.e., Cee — 0.09172(1 + 72). (As the density increases, 
the maximum in p{T) becomes less pronounced, and 
the non-monotonicity eventually disappears at rtg > 
1 X 10^^ cm^^, as the ballistic region is approached.) 
The solid lines are the universal theoretical curves 
for p(r)//9max and 72(7) plotted as a function of 
Pi„axln(T/Tmax)- One sees from Fig.Hthat at T = Tmax, 
the value of the interaction amplitude 72 ~ 0.45, which 
is in remarkable agreement with theory for n„ = 2. The 
agreement between theory and experiment is especially 
striking given that the theory has no adjustable param- 
eters. Systematic deviations from the universal curve 
occur at lower densities as higher order corrections in p 
become important. 

The data for 72 allow us to calculate the renormalized 
Lande g- factor g* — 2(1-1-72). For our highest density, it 
grows from g* « 2.9 at the highest temperature to g* « 4 
at the lowest. This is the first time the g-factor has been 
measured in the diffusive regime and is seen to increase 
with decreasing temperature. Earlier experiments were 
done in the ballistic regime where the g-factor was found 
to be g* « 2.8 and temperature independent [isj . 

Finally, we would like to note that the 2D electron sys- 
tem in silicon, used in this study, has a short-range disor- 
der potential and constitutes a multi-valley system: fea- 
tures assumed in the theory. In other material systems, 
including GaAs/AlGaAs heterostructures and SiGe het- 
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FIG. 4; Comparison between theory (lines) and experiment (symbols), (a): p//5max as a function of pmax ln(r/rniax). (b): 
72 as a function of pmax ln(r/rniax). Vertical dashed lines correspond to T = Tmax, the temperature at which p{T) reaches 
maximum. Note that at this temperature, the interaction amplitude 72 « 0.45 (indicated by the horizontal dashed line in 
(b)), in excellent agreement with theory. Symbols are the same as in Figs [T] and |3] and correspond to Us — 9.87, 9.58 and 
9.14 X 10^° cm-^ 



erostructures, the disorder potential is long-range, which 
leads to very high mobilities. The diffusive regime in 
these systems is therefore hard to reach at reasonable 
temperatures. In addition, high mobilities imply very low 
carrier densities, which leads to a plethora of finite tem- 
perature corrections arising from non-degeneracy effects, 
large-r^ screening effects, etc. Therefore, the transition 
observed so far in these systems (see, e.g., Refs. 3, 
might be a "high" -temperature crossover phenomenon. 
Ultimately, however, at low enough temperatures the dif- 
fusive regime will unavoidably be reached. It remains to 
be seen whether or not the low temperature behaviour 
in other 2D systems is quantitatively described by the 
two-parameter scaling theory. 
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